A new mechanism of hypersensitive transport in tilted sharp ratchets 

induced by noise flatness 
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The noise-flatness-induced hypersensitive transport of overdamped Brownian particles in a tilted 
ratchet system driven by multiplicative nonequilibrium three-level Markovian noise and additive 
white noise is considered. At low temperatures the enhancement of current is very sensitive to the 
applied small static tilting force. It is established that the enhancement of mobility depends non- 
monotonically on the parameters (flatness, correlation time) of multiplicative noise. The optimal 
values of noise parameters maximizing the mobility are found. 



PACS number(s): 05.40.-a, 05.60.Cd, 02.50.-r 

Recently, noise- induced hypersensitivity to small time- 
dependent signals in nonlinear systems with multiplica- 
tive noise has been the topic of a number of physical 
investigations [1-4]. A motivation in this field has come 
from numerical, analytical, and experimental studies of 
a nonlinear Kramers oscillator with multiplicative white 
noise [1,4]. Under the effect of intense multiplicative 
noise, the system is able to amplify an ultrasmall de- 
terministic ac signal (of the order of, e.g., 10~ 20 ) up to 
the value of the order of unity [1]. Afterwards, a re- 
lated phenomenon such as noise-induced hypersensitive 
transport was found in some other systems with multi- 
plicative dichotomous noise [2,3]. Noise-induced hyper- 
sensitive transport was also established in a phase model, 
i.e., dep/dt — a — b sirup, with a strong symmetric mul- 
tiplicative colored noise. It was shown that in such a 
system a macroscopic flux (current) of matter appears 
under the effect of ultrasmall dc driving [2] . It is impor- 
tant to notice that the physical mechanism underlying 
the phenomenon of hypersensitive transport presented in 
Refs. [2,3] is based on the assumption that the periodic 
potential is smooth. It is easy to see that in the case of 
a periodic sharp potential the above mechanism cannot 
bring forth any hypersensitive transport. 

Theoretical investigations [5-9] indicate that noise- 
induced nonequilibrium effects are sensitive to noise flat- 
ness, which is defined as the ratio of the fourth moment 
to the square of the second moment of the noise pro- 
cess. Although its significance is obvious, the role of the 
flatness of fluctuations has not been researched to any 
significant degree to-date. In the present paper we as- 
sume the multiplicative noise to be a zero-mean trichoto- 
mous Markovian stochastic process [10]. It is remark- 
able that for trichotomous noises the flatness parameter 
ip, contrary to the cases of the Gaussian colored noise 
(ip = 3) and symmetric dichotomous noise (ip — 1), can 
have any value from 1 to oo. The flatness as an ex- 
tra degree of freedom (in comparison with dichotomous 
noise) can prove useful when modelling actual fluctua- 
tions, e.g., thermal transitions between three configura- 
tions or states. This is the reason why we choose in 
the phase space of possible non-equilibrium models the 



trichotomous noise. Although both dichotomous and tri- 
chotomous processes may be too rough approximations 
of the actual non-equilibrium fluctuations, the latter is 
more flexible, including all cases of dichotomous pro- 
cesses and, as such, revealing the essence of its peculiar- 
ities. A further virtue of the models with trichotomous 
noise is that they constitute a case admitting exact ana- 
lytical solutions for some nonlinear stochastic problems, 
such as colored-noise-induced transitions [10] and rever- 
sals of noise-induced flow [9] . 

The main purpose of this paper is to establish a new 
mechanism of hypersensitive transport, demonstrating 
that the flatness of multiplicative noise can generate hy- 
persensitive response to the small external static force 
in a tilted sharp ratchet system. We will show that in 
the region of hypersensitive response the value of mobil- 
ity can be controlled by means of thermal noise. For low 
temperatures, we find that the mobility exhibits resonant 
behavior at intermediate values of the parameters of the 
multiplicative noise (flatness, correlation time). 

We consider an overdamped multinoise tilted ratchet, 
where particles move in a one-dimensional spatially pe- 
riodic potential of the form V(x,t) = V(x)Z(t), where 
Z(t) is a trichotomous process [10] and V(x) is a piece- 
wise linear function, which has one maximum per period. 
The additional force consists of thermal noise with tem- 
perature D, and an external static force F. The system 
is described by the dimcnsionlcss Langevin equation 
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Z(t)h{X)+F + i{t), h(x) 



dV(x) 
dx 



(1) 



where V(x) = V(x)/Vo, V(x) is a spatially periodic func- 
tion with period L and Vq = V rnax — V m in ■ The usual di- 
mensionalized physical variables are indicated by tildes 
and the space and time coordinates read X = X/L 
and t — tVo/nL 2 with k being the friction coefficient; 
F = VqF/L is a constant external force. The thermal 
noise satisfies (£(*)) = and (£(*i)£(t 2 )) = 2D5(t 1 - t 2 ). 
Regarding the random function Z(t), we assume it to 
be a zero-mean trichotomous Markovian stochastic pro- 
cess [10] which consists of jumps among three values 
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z = {1,0, —1}. The jumps follow in time according to a 
Poisson process, while the values occur with the station- 
ary probabilities = P s (—1) = q and P s (0) = l — 2q. 
In a stationary state, the fluctuation process statisfies 
(Z(t)) = and (Z(t + r)Z(t)) = 2qcxp(-vr), where the 
switching rate v is the reciprocal of the noise correlation 
time t c = \jv. The trichotomous process is a special 
case of the kangaroo process [5] with a flatness parame- 
ter <p = (Z 4 (t))/(Z(t)} 2 = l/(2q). At large flatnesses our 
trichotomous noise essentially coincides with the three- 
level noise used by Bier [6] and Elston and Doering [7] . 
The master equation corresponding to Eq. (1) reads 

^ t Pn(x,t) = - — [r n P n (x i t)]+J2UnrnP m (x,t) , (2) 

m 

where T„ = z n h(x)+F — Dd x and P n (x,t) is the 
probability density for the combined process (x, z n ,t), 
n,m — 1,2,3, Z\ = 1, z% = 0, z 3 = — 1, and 
Uik = v\q + (1 — 3g)(5i2 — &ik\- The stationary current 
J = ^2 n j n (x) is then evaluated via the current densi- 
ties j n (x) = (z n h(x) + F - Dd x )P^{x), where is 
the stationary probability density in the state (x,z n ). 
To calculate the stationary probability density in the 
x space, P{x) = J2 n Pn( x ), anc > the stationary cur- 
rent, J — const, six conditions are imposed on the so- 
lutions of Eq. (2), namely, the conditions of periodic- 
ity Pn(x) = Pn( x + 1)) n = 172,3, and normalization of 
P£ (x) over the period interval L = 1 of the ratchet poten- 
tial Z(t)V(x), which read J Q P({x)dx — L P^(x)dx = q 

and P$(x)dx = 1 - 2q. 

To derive an exact formula for J, we assume that the 
potential Z(t)V(x) = Z{t)V(x - 1) in Eq. (1) is piece- 
wise linear (sawtoothlike) and its asymmetry is deter- 
mined by a parameter d £ (0, 1), with V{x) being sym- 
metric when d = 1/2. A schematic representation of 
the three configurations assumed by the " net potentials" 
V n (x) = z n V(x) — Fx associated with the right hand side 
of Eq. (1) is shown in Fig. 1. Regarding the symmetry of 
the dynamic system (1), we notice that J(—F) = —J(F) 
and J(F, d) = J(F, 1 — d). Thus we may confine ourselves 
to the case d < 1/2 and F > 0. Obviously, for F = 0, 
the system is effectively isotropic and no current can oc- 
cur. In the case of zero temperature the both noise levels 
z n=1 , 3 = ±1 in Eq. (1), where F < min{l/d, 1/(1 - d)}, 
give zero flux. However, if one allows switching between 
the three dynamic laws V n (x), n = 1,2,3, the resulting 
motion will have a net flux which can be much greater 
that the flux by the dynamic law V2 = —Fx. If the 
rate of reaching the minimal energy in each well consid- 
erably exceeds the switching rate v, the leading part of 
the net flux is achieved in the following way: a parti- 
cle locked in the potential minimum 1 switches to point 
2, then slowly moves to point 3, switches to point 4 (or 
to 5 with equal probability), and rapidly slides down to 
point 6 (or from 5 back to 1), etc (see Fig. 1 and cf Ref. 
[11]). In this case hypersensitive transport is possible and 
can be intuitively understood. The described scheme is 
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FIG. 1. The mechanism of hypersensitive transport. The 
lines depict the net potentials V„(x) — z n V(x) — Fx with 
zi = 1, Z2 = 0, and 23 = — 1. A particle cannot move of its 
own accord along the potentials V\ and V3. However, if one 
allows switching between the potentials V n , n — 1,2,3, the 
particle will move downhill along the trajectory 1:2:3:4:6. 

valid only in the absence of thermal noise. Otherwise, a 
particle is able to pass by a thermally activated escape 
across the potential barriers in both directions. However, 
it predominantly moves to the right and hypersensitive 
transport still occurs (at least at sufficiently low tem- 
peratures). As the "force" h(x) = —dV(x)/dx is piece- 
wisely constant, h(x) = h\ = 1/d for x £ (0, d) (modi) 
and h(x) = h 2 = -1/(1 - d) for x £ (d, 1) (modi), 
Eq. (2) splits up into two linear differential equations 
with constant coefficients for the two vector functions 
P?(x) = P&) (i = 1,2) defined on the intervals 

(0, d) and (d, 1), respectively. The solution reads 

5 

PZi(x) = JA n + Y,C lk A mk e x ^ x , (3) 
fc=i 

where Cik are constants of integration, the constants A n 
and A nik are given by A x = A 3 = qJ/F, A 2 = (1 - 
2q)J/F, A mk = (DX ik - F)[D\% — (F — z n hi)X ik - v] 
for n = 1,3, A 2ik = 2/if X ik - (Auk + A 3ik ), and {X ik ,k = 
1, . . . , 5} is the set of roots of the algebraic equation 

D 3 Xl - 3D 2 FXf + D(3F 2 - 2Dv - ftf)Af 
+ F(ADv - F 2 + h 2 )X 2 + v{Dv - 2F 2 + 2qh 2 )X l 

- v 2 F = . (4) 

Eleven conditions for the ten constants of integration of 
Eq. (3) and for the probability current J can be deter- 
mined at the points of discontinuity, by requiring con- 
tinuity, periodicity, and normalization of P?(x). This 
procedure leads to an inhomogeneous set of eleven lin- 
ear algebraic equations. Now, an exact formula for the 
current J can be obtained as a quotient of two determi- 
nants of the eleventh degree. The exact formula, being 
complex and cumbersome, will not be presented here, 
however, it will be used to find (i) the dependence of the 
current J on the tilting force F and the dependence of 
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the mobility ra — J/F on the flatness if — l/(2q), which 
are displayed in Figs. 2 and 3, respectively, and (ii) the 
asymptotic limits of the current J at low temperature 
and small external force. 

Figure 2 shows the induced current J as a function of 
the external force F for two different values of tempera- 
ture and for three different values of d with fixed ip = 2.5 
and v = 8. In this figure, one also observes the hypersen- 
sitive response at very low forcing, which apparently gets 
more and more pronounced as the thermal noise strength 
D decreases. For the case D = 0, d = 0.5, the results of 
Monte Carlo simulations of the current J = J(F) are 
also presented. The tendency apparent in Fig. 2, namely 
a decrease in the mobility for very low forcing as the 
asymmetry of the potential grows, is also valid for large 
asymmetries, e.g., when d < 0.05. 

To obtain more insight, we shall now study some 
asymptotic limits of the current. 

At the fast-noise limit, we allow v to become large, 
holding all the other parameters fixed. Thus, at very high 
frequencies of colored fluctuations, the system is under 
the influence of the average fluctuating potential. In the 
v — > oo limit, the current is then given by 

J = F + 0{v~%) . 

The form of the leading term of the current J is not 
confined to the fast-noise limit. It is also valid for the 
asymptotic limit of a high temperature, D » 1, and in 
the case of a large "load" force F (F — > oo, all the other 
parameters fixed). 

At the long-correlation-time limit v — > 0, the equations 
(2) for -Pf(x), Pi{x), and Pi(x) are decoupled and the 
total current is given by the average of each current for 
the corresponding potential configurations. In the case of 
the symmetric potential, d = 1/2, the current J saturates 
at the value 
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2q)F 
2q(4 



F 2 ) 2 sh (F/2D) 



16D[ ch i - ch - F(4 - F 2 ) sh 
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For F < 2, we can see that the current J tends to 
(1 — 2q)F as D — > 0. This result is consistent with the 
physical intuition that the probability densities Pi(x) 
and P%(x) are ^-distributed at deterministic stationary 
states (minimums of potentials): the random variable Z 
takes values ±1 for a sufficiently long time to allow the 
deterministic stationary state to be formed. 

In the case of zero temperature, D — 0, and symmetric 
potential, d = 1/2, one finds from the exact formula that 
on the assumption F < 2 the current equals 



J = vF 



A X C 2 - C X A 2 
B\C 2 — B 2 Ci 



(5) 



where A, = F{a t [F - (4 - F 2 ) Vl } - 2(1 - 2q)}, B t = 
(v + 16q)A l + 32 9 (1 - 2q)(2a. l + F), C t = qA t + 2(1 - 
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FIG. 2. The current J vs applied force F in the region of 
the hypersensitive response. The flatness parameter equals 
tp = 2.5 and the switching rate v = 8. Solid straight line: 
D = 0, d = 0.5. Dotted line: D = 4 x 10" 8 , d = 0.5. Dashed 
line: D = 4 x 10" 8 , d = 0.2. Solid curved line: D = 4 x 10" 8 , 
d — 0.05. The filled dots on the solid straight line are obtained 
by means of Monte Carlo simulations. Notice the jump of the 
current from the zero level to the solid line corresponding to 
the infinite derivative of J(F) at F — 0. 

2q)[4r li + Fq + 2a i (l + Fr] i )}, rji = F _1 (4 — F 2 )~ x [F 2 - 
Aq - 2e t (4q 2 + F 2 (l - 2q)) 1 / 2 }, a t = th(^/4), * = 1,2, 
£1 = 1, £2 = -1. 

Thus, at the low-force limit, F — > 0, the current will 
saturate at the finite value 



lim J 

F->0 



Ja 



32iyq(l - 2q) 
(*/ + 8) 2 



(6) 



As J(F = 0) = 0, the hypersensitive response is ex- 
tremely pronounced in this case, with the current picking 
up with an infinite derivative at F — (see also Fig. 2). 
The asymptotic current J a exhibits a bell-shaped (reso- 
nance) form as v or q are varied. The optimal correlation 
time T m that maximizes the current equals 1/8, and the 
optimal flatness ip m — l/(2q m ) = 2. It is remarkable that 
in the case of a dichotomous noise q = 1/2, the hyper- 
sensitive response disappears and in the low forcing limit 
the leading-order term of the current is proportional to 
F : J re vF(u+ 12)/(z/ + 8) 2 . 

At the low forcing limit, F <C 1, a natural way to 
investigate the behaviour of J is to apply small- -F per- 
turbation expansions. A stationary solution of Eqs. (3) 
and (4) with D 7^ 0, d— 1/2 is constructed in terms of 
integer powers of F. The current can be expressed as 
,/ = Fnii + F 2 ni2 + ■ ■ ■■ We shall calculate the leading 
term of the current Fm\. Notably, the analysis of this 
section is valid for the values of parameters satisfying the 
condition F < (2qvD) 1 ^ 2 . This condition results from 
the assumption that the higher-order terms in the expan- 
sion of the roots of Eq. (4) are asymptotically smaller 
than the lower-order terms held in the calculation. At 
sufficiently small temperature, D <C min{l, 2qv 1 8q/i>}, 
the formula for the leading-order term Fm\ of the cur- 
rent is 
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FIG. 3. The mobility m = J/F vs the flatness pa- 
rameter q = l/(2<p) at d = 1/2, D = 4 x 10~ 8 , 
and F = 10~ J . The curves computed from the ex- 
act formula for the current J correspond to the values of 
v = 8/3, v = 1, z/ = 8, v = 0.1, ^ = 100 from top to bottom. 
The non-monotonic sequence of the values of v stems from the 
bell-shaped dependence J = J(y). Note that the maximum 
of the mobility lies at q — 1/6 and v — 8/3. The dots were 
computed by means of the asymptotic formula (7). 
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Here the symbol G stands for the terms which do not 
increase as D — > 0. An extreme sensitivity of the mobility 
m to thermal noise can be seen from the factor D^ 1 ^ 2 in 
Eq. (7) that increases unboundedly as D — > 0. It can be 
seen easily that the functional dependence of the mobility 
on the flatness ip and on the correlation time r c is of a 
bell-shaped form. The mobility mi reaches a maximum 
at the flatness ip m — 3 and at the correlation time r m = 
3/8. The dependence of the mobility m = J/F on the 
parameters q and v for a fixed force value F = 10~ 5 and 
for a fixed temperature D = 4x 10~ 8 is shown in Fig. 3. 
We can see that the asymptotic formula (7) is in excellent 
agreement with the exact results. 

Let us note that the sufficient condition, F < ^J2qvD, 
has a distinct physical mea ning: the characteristic dis- 
tance of thermal diffusion y D/v is larger than the typi- 
cal distance F/v for the particle driven by the determin- 
istic force F in the state z — of the trichotomous noise. 
Let us look at the latter statement more closely on the as- 
sumption that v -C 1. For this assumption within the in- 
terval (0, 1), the probability distributions P^(x), n = 1,3, 
are, evidently, concentrated at x — (or x — 1/2). Next, 
we shall consider the trajectory (1 : 2 : 3 : 4 : 6) in 
Fig. 1. The particles locked at the potential minimum 1 
(x = d = 1/2) will go at the initial time t = to point 
2, where z = 0. The first time when the noise turns to 
either z = 1 or z = — 1 is denoted by to- As the time 
of movement from 4 to 6 is much less than to, it is easy 
to find that during the time interval (0, to) the center of 
mass has shifted by 
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2^/^Dt 



I rFto 



exp< - 




In the case of a trichotomous noise the probability 
W(t) that in a certain time interval {0,t) the tran- 
sitions z = 0^z = ±1 do not occur, is given by 
W(t) = exp(— 2qvt). The probability that the transition 
z = 0^> z = — 1 occurs within the time interval (t, t + dt) 
is qvdt. Consequently, 



(Ax) =qv I c- 2qi,t0 A x dto 
Jo 



F 



8\/2^7D 



Considering that the average number of transitions per 
unit of time into the state z = is 2qv(l — 2q), we ob- 
tain J = 2qv(l - 2<?)(Ax)«F(l - 2q)^qDj%\fD. Thus, 
we have obtained an earlier result, namely, Eq. (7) for 
v <C 1. Formula (7) is one of our main results. Note 
that the above procedure can be repeated in a straight- 
forward but tedious way for more complicated cases in- 
volving asymmetric potentials and potentials with several 
extrema per period. The phenomen is robust enough to 
survive a modification of the multiplicative noise. The 
key-factor is the noise flatness, indicating how long the 
noise level dwells on the state z = 0. If the flatness 
parameter is greater than 1, the effect does exist. For 
example, the multiplicative noise can also be a Gaussian 
stationary process. 

It is quite remarkable that the above results are also 
applicable for amplifying adiabatic time-dependent sig- 
nals, i.e., signals of much greater periods than the char- 
acteristic time of establishing a stationary distribution, 
even in the case of a small input signal-to-noise ratio 

f/Vd < i. 

We emphasize that our mechanism of hypersensitive 
transport is of a qualitatively different nature from a 
recently found effect, where a noise-induced enhance- 
ment of the current of Brownian particles in a tilted 
ratchet system has also been established [2,3]. In the lat- 
ter scenario, a system with a periodic smooth potential 
exhibits hypersensitivity under the effect of multiplica- 
tive dichotomous noise because of noise-induced escape 
through fixed points of the dynamics. This occurs be- 
cause the stable and unstable fixed points of the alterna- 
tive dynamics, which coincide in the absence of the tilt 
F, are shifted apart by a small force (see also Ref. [11]). 
In the mechanism reported here, we have a sharp peri- 
odic potential: the stable and unstable fixed points of the 
dynamics coincide also for any small tilt. The crossing of 
the location of the fixed points is achieved by a combined 
influence of the flatness of the multiplicative noise and a 
small tilt-forcing. 

In a general case, if the potential is smooth and the flat- 
ness of multiplicative noise is greater than 1 , both mech- 
anisms play an important role and should be taken into 
account. Our calculations show that the factor Fy/ v/D 
in Eq. (7) is generated by thermal diffusion in the state 
= 0, while the circumstance that the potential is sharp 
has no effect on this factor. On the other hand, for adi- 
abatic switching, the mechanism described in Ref. [2] 
generates the current J <~ vFj\f~D. Consequently, our 
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mechanism for sufficiently small switching rates induces [12] 0. V. Gerashchenko, S. L. Ginzburg, and M. A. Pus- 
hypersensitive transport more effectively than the one tovoit, Eur. Phys. J. B 19, 101 (2001). 
proposed by Ginzburg and Pustovoit. This conclusion is 
in agreement with the results of Ref. [2] , presenting nu- 
merical simulations of the phenomenon of hypersensitive 
transport based on a phase model with the multiplica- 
tive colored Gaussian noise (ip — 3). It is established 
that in the case of low switching rates the transport for 
the Gaussian noise appears to be more effective than for 
dichotomous stimuli. Regrettably, the authors of Ref. [2] 
did not consider the role of noise flatness and the physics 
of this discrepancy. 

In conclusion, the reported mechanism of generating 
hypersensitive transport by the flatness of multiplicative 
noise is of general relevance for many physical, biological, 
and chemical systems, and may provide another possi- 
bility to control signal amplification. The sensitivity of 
system response to small input signals can be either en- 
hanced or suppressed by changing the noise parameters 
(flatness, correlation time, temperature). In agreement 
with Ref. [12], we believe that the phenomenon proposed 
may also shed some light on the ability of biological sys- 
tems to detect weak signals in a noisy environment. 

We acknowledge partial support by the Estonian Sci- 
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